Local dynamical lattice instabilities: prerequisites for resonant pairing 

superconductivity 



00 
O 

o 

00 



o 
o 

3 



I 

C 

O 

o 



> 
in 
in 

ON 

o 

oo 
o 



- 1—^ 

X 



Julius Ranninger 1 and Alfonso Romano 2 
1 Institut Ned, CNRS-UJF, BP 166, 38042 Grenoble, France 
2 Laboratorio Regionale SuperMat, CNR-INFM, Baronissi (Salerno), Italy, and 
Dipartimento di Fisica "E. R. Caianiello" , Universita di Salerno, 1-84081 Baronissi (Salerno), Italy 

(Dated: September 8, 2008) 

Fluctuating local diamagnetic pairs of electrons, embedded in a Fermi sea, are candidates for non- 
phonon-mediated superconductors without the stringent conditions on T c which arise in phonon- 
mediated BCS classical low-T c superconductors. The local accumulations of charge, from which such 
diamagnetic fluctuations originate, are irrevocably coupled to local dynamical lattice instabilities and 
form composite charge-lattice excitations of the system. For a superconducting phase to be realized, 
such excitations must be itinerant spatially phase-coherent modes. This can be achieved by resonant 
pair tunneling in and out of polaronic cation-ligand sites. Materials in which superconductivity 
driven by such local lattice instability can be expected, have a T c which is controlled by the phase 
stiffness rather than the amplitude of the diamagnetic pair fluctuations. Above T c , a pseudogap 
phase will be maintained up to a T* , where this pairing amplitude disappears. We discuss the 
characteristic local charge and lattice properties which characterize this pseudogap phase and which 
form the prerequisites for establishing a phase-coherent macroscopic superconducting state. 

PACS numbers: 71.38.-k, 74.20.Mn, 72.20.Jv 



I. INTRODUCTION 

Ever since the BCS theory identified the phonon ex- 
change mechanism as the cause for electron pairing lead- 
ing to a superconducting state and, at the same time, 
pinned down the order of magnitude for the critical tem- 
perature T c , considerable efforts were undertaken to by- 
pass these constraints on T c . One of the leading ideas 
was to look for fluctuating local diamagnetism caused by 
intrinsic atomic correlations which would act as an es- 
sential component for pairing in the many-body electron 
wavefunction. Since diamagnetic pairs of electrons nec- 
essarily lead to local lattice deformations, one searched 
for materials* where such deformations would exist on 
a dynamical level, as close as possible to, but not at, a 
global lattice instability. This was thought to favor op- 
timal diamagnetism and strong superconducting correla- 
tions, without leading to localization which would origi- 
nate from static lattice instabilities. 

In this line of thinking, systematic studies were un- 
dertaken searching for compounds liable to show charge 
disproportionation, originating from unstable "formal" 
valence states of the cations which build up such struc- 
tures. Those formal valence states, required by the 
chemical stoichiometry given a certain crystalline envi- 
ronment, are skipped. As a result, such materials be- 
come composed of a mixture of cations with one more 
and one less charge, respectively. Examples for such va- 
lence skippers which are unstable are Tl 2+ , Pb 3+ , Sn 3+ , 
Bi 4+ , Sb 4+ . These elements prefer to exist as Tl 1+ , 
Pb 2 +, Sn 2 +, Bi 3 +, Sb 3 + together with Tl 3 +, Pb 4 +, Sn 4 +, 
Bi 5+ , Sb 5+ . It implies a tendency to electron pairing 
and for that reason has attracted, early on, much at- 
tention for possible non-phonon-mediated superconduc- 
tivity of purely atomic origin. Charge disproportion- 



ation then originates from effective negative-C/ centers 
to which the electrons, associated with the cations with 
unstable valence configurations, are attracted and form 
bound electron pairs. Casting such physics in terms 
of a negative-?/ Anderson model gives rise to a vari- 
ety of macroscopic phases, such as charge fluctuation- 
driven superconductivity^!^ correlation-driven insu- 
lating states, translationally symmetry-broken charge- 
ordered states, etc. For the exactly half-filled band case, 
the stable valence configurations (differing by two elec- 
trons) are degenerate and as a result ensure resonant pair 
tunneling in and out of charge Kondo impurities. 

The mechanism behind valence skipping has been 
a highly debated subject and for a long time was 
attributed^' 7 -'^ to an intrinsic intra-atomic interplay be- 
tween the electron affinity and the ionization energy, 
combined with a crossing of stable closed 4s, 5p, 6s elec- 
tron shells. It has become increasingly evident in recent 
years that valence skipping and the resulting charge dis- 
proportionation systematically appear together with dis- 
tinct deformations of the cation-ligand complex, which 
strongly depend on which valence state they are in. 
Thus, in order to overcome the relatively large on-site 
Coulomb repulsion, 1 ^ strong Madelung potentials and 
strong covalent bonding ) 11 : 12 an adequate polarizibility 
of the material and sufficiently large dielectric constants 
(~ 20 — 40 ^ must play together so that charge dispro- 
portionation can be realized. It is now largely consid- 
ered that the prime mechanism for charge disproportion- 
ation lies in the local lattice relaxations of the ligands 
surrounding those cationsJ^ This picture, closer to the 
traditional scenario originally hypothesized for localized 
defects in semiconductors,— however now requires tak- 
ing into account a time retarded - rather than static - 
exchange interaction between bound electron pairs and 
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itinerant electrons. 

The question we address in this paper is to what ex- 
tent charge disproportionation can be realized in form of 
dynamical diamagnetic double charge fluctuations on ef- 
fective sites, composed of cations and their surrounding 
ligand environments. The latter are deformed by local 
lattice instabilities which change their bond-lengths or 
angles. In this way, two electrons are either captured 
from the Fermi sea in form of a self-trapped bipolaron or 
are released onto the cations in the immediate vicinity 
where they become itinerant again. One can thus hope 
for a resonant tunneling of pairs of electrons in and out of 
such polaronic sites, inducing local dynamical pair cor- 
relations in the itinerant sector of the electron system. 
In this sense our scenario is similar to that of fermionic 
atomic gases in which a global supcrfiuid phase is induced 
by a Feshbach resonance^ which describes the charge ex- 
change interaction between bound electron-spin triplet 
states and electron-spin singlet scattering states. For 
such coherent pair tunneling to occur in our counterpart 
solid state example, the charge and ligand deformation 
fluctuations will have to be highly correlated with each 
other. We shall examine here to what extent this is fea- 
sible. If the exchange of the two types of double charge 
carriers is sufficiently efficient, such charge disproportion- 
ation materials and, possibly in a wider sense, the copper 
oxides and their non-BCS like superconductivity, could 
be due to the resonant bipolaron mechanism described 
above. 

Particularly well studied examples of such materials 
are: 

(A) Pbi-zTlzTeiS with the parent compound PbTe 
which is a small gap semiconductor. Upon doping with 
Tl, which exist as Tl 3+ in this compound and plays the 
role of a kind of negative-?/ center, Pbi_ x Tl x Te becomes 
a superconductor with a relatively high value of T c ~ 
1.5 K, considering its small carrier density. For details we 
refer the reader to an extensive literature on this subject, 
see refs. 7, 8 and references therein. 

(B) BaBi^Pbi-^Oa— and Beti-^K^BiOa with the par- 
ent compound BaBi03 which is a diamagnetic insulator 
with a charge-ordered state involving alternating Bi m 
and Bi v (from now on we shall use the more appropri- 
ate denomination accounting for covalency for the formal 
valence states, indicated by roman numeral in the super- 
scripts). Band theory predicts a metal on the basis of 
a formal valence state Bi lv . The charge-ordered state 
is stabilized due to a highly anisotropic polarization of 
the O ligand environment^ inciting positional changes 
in the two stable valence species of the Bi ions. Bi v is 
in a regular octahedra ligand environment, with a Bi— O 
distance of 2.12 A. Bi 111 is in a pseudo-octahedral lig- 
and environment where one of the six oxygen ions in the 
octahedra is displaced to such an extent (with a corre- 
sponding Bi— O distance of 2.28 A) that it is effectively 
becoming 2 ~ after having transferred an electron to the 
cationJ^ The total outcome is a structurally different lig- 
and environment which favours a Bi 111 valence state. Par- 



tially substituting Ba by K results in a superconducting 
state with a relatively high value of T c ~ 30 K, which 
has been considered as an example of superconductivity 
induced by double-valence fluctuations^^ 

(C) Charge disproportionation is also seen on a molec- 
ular, rather than atomic, level in systems such as 
T1407) 1 - V 2 0g>22> and W03-K- 1 in which one observes 
diatomic pairs of Ti m -Ti m together with Ti 4 +-Ti 4+ , 
V IV -V IV together with V v -V v , and W v -W v to- 
gether with W VI — W VI , respectively. These different 
valence states of the molecular units again show dis- 
tinctly different intra-molecular distances, such as 2.8 A 
for Ti m -Ti m and 3.08 A for Ti IV -Ti IV bonds. More- 
over, with increasing temperature these systems exhibit 
a thermally activated conductivity by pair hopping. This 
demonstrates the stability of such bound electron pairs 
on deformed molecular units which are locked together 
in a dynamical fashion with the charge fluctuations. 

(D) Indications for dynamical charge fluctuations are 
also seen in the cuprate superconductors. Here the situ- 
ation is however more involved. Stable valence states of 
Cu exist as: 

(i) Cu 1 in a dumbbell oxygen ligand environment, with 
a characteristic Cu— O bond-length of 1.84 A (examples 
are Cu 2 0, BaCu 2 2 , FeCu0 2 ); 

(ii) Cu 11 in square planar oxygen ligand environments 
with a characteristic Cu— O bond-length of 1.94 A (an 
example is CUO4); 

(iii) Cu 111 in square planar oxygen ligand environ- 
ments, in rare cases, with a characteristic Cu— O bond- 
length of around 1.84 A, and hence practically identical 
to that of Cu 1 in dumbbell environments 1 ^ (an example 
is KCu0 2 ). 

In the cuprates the formal valency of Cu lies between 
II and III. The square lattice of the cuprates is neither a 
strict square planar nor a dumbbell ligand environment 
and there is no indication for a static disproportionation 
of Cu 11 . Cu 111 , which should be favored in this cuprate 
planar lattice structure, could appear only through rapid 
metallic fluctuations or be due to the stoichiometry of the 
structure, with an atypical Madelung potential. 22 Yet, 
local dynamical lattice fluctuations involving the oxygen 
ligands of the planar Cu ions, the Cu— O bond stretch 
(or buckling) modes, exist. They can be considered as 
being linked to charge fluctuations involving local ligand 
environments which are typical of stable Cu 1 and Cu 111 . 
In the anti-nodal region of the Brillouin zone, where the 
electronic density of states (DOS) is dominated by the the 
pseudogap phenomenon, local scanning tunneling micro- 
scope (STM) spectroscopy 23 did establish a strong cou- 
pling of electrons with wave vectors [k Xl k y ] = [±7T, 0] and 
[0, ±7r] and local phonon modes with characteristic fre- 
quencies distributed around 50 meV. The isotope varia- 
tion of these modes (most likely associated to the Cu— O 
bond stretch mode) correlates with that of the pairing 
gap. This strongly indicates that this local Cu— O bond 
stretch mode could play a significant role in the pairing 
mechanism of the cuprates. 
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If double charge fluctuations do play a significant role 
in establishing a superconducting phase, they have to 
sustain a free particle- like behavior with spatial phase co- 
herence. This should occur in spite of the sizeable lattice 
relaxation which, a priori, tends to render them rather 
diffusive. If in this largely diffusive dynamics we can 
obtain a finite component of coherent phase-correlated 
double charge fluctuations, then a superconducting state 
can materialize. Such a state does not describe a purely 
bosonic system of bipolarons and bipolaronic supercon- 
ductivity, but is characterized by resonating bipolarons 
embedded in a Fermi sea. Their signature is a pseudo- 
gap in the single-particle DOS at the Fermi level due 
to strong pairing correlations above the superconducting 
phase. This scenario has been discussed in great detail 
over the past years on the basis of the Boson-Fermion 
Model by J. Ranninger and collaborators. 

To address the question of the feasibility of coherent 
pair tunneling, we consider charge carriers capable of 
existing either as selftrapped bipolarons or as itinerant 
quasi-free tight-binding charge carriers on a small clus- 
ters consisting of: 

(i) cation-ligand complexes which are deformable and 
which can capture the charge carriers in form of bipo- 
larons; 

(ii) square plaquettes composed of four structurally 
identical cations-ligand complexes (surrounding the de- 
formable central cation-ligand complex), which we con- 
sider as undeformable. Charge carriers on these plaque- 
ttes shall hence move as itinerant entities rather than 
being captured by polaronic effects. 

Double charge fluctuations in the present picture im- 
plies fluctuations between doubly occupied and unoccu- 
pied sites, driven by charge carriers tunneling in and out 
of such sites and ultimately resulting in pair correlations 
among those on the plaquette sites. Such essentially lo- 
cal physics contains the essence of triggering a crossover 
from an insulating purely charge-disproportionated state 
to a superconducting state upon tuning certain parame- 
ters such as the electron-lattice coupling, the adiabaticity 
ratio or the charge carrier density. This physics, intrinsi- 
cally related to a very local mechanism^! manifests itself 
in the single particle spectral function, the local diamag- 
netism, the local phonon softening, the quasi-elastic peak 
in the neutron scattering cross-section for the phonons 
and the double-peaked pair distribution function for the 
bond-length fluctuations. 

In section II we discuss the salient features of such 
double charge fluctuations driven by bond-length fluctu- 
ations and present a scenario in terms of a kind of Hol- 
stein model which takes into account the different ligand 
fluctuations characterizing the various charge configura- 
tions. In section III we propose a model Hamiltonian 
which we consider to be adequate to describe the local 
physics of such systems. We discuss its basic features as 
far as the efficiency of local pair tunneling is concerned 
and how such dynamically correlated charge-ligand de- 
formation fluctuations are locked together in order to 



form phase-coherent correlated itinerant excitations on a 
short spatial scale. In section IV we discuss the specific 
characteristics of: (i) the local lattice dynamics close to 
instabilities, phonon softening, the quasi-elastic peak in 
the dynamical structure factor and the double peak struc- 
tured pair distribution function, (ii) the diamagnetic cor- 
relations induced among the itinerant electrons via such 
dynamical local lattice instabilities and how they mani- 
fest themselves in a pseudogap feature in the electronic 
density of states. In the Summary, section V, we briefly 
review the main findings of the present work and discuss 
the further issues which have to be treated next. 



II. ELECTRON LOCALIZATION VERSUS 
DELOCALIZATION IN ELECTRON-LATTICE 
COUPLED SYSTEMS 

The crossover from metallic, or superconducting, to lo- 
calized behaviour in electron-lattice coupled systems has 
been a topic of longstanding interest and debate.— The 
issue has been reinvestigated with renewed vigor in re- 
cent years, when it became apparent that an insulating 
phase of localized polaronic charge carriers, bordering on 
a superfluid phase, can be induced by many-body cor- 
relation effects . 26 ' 27 ! 28 Earlier widely pursued ideas con- 
sidered density-driven stripping off of the phonon clouds 
of individual polaron o 29 ' 30 and a related to it possible 
breakdown of Wigner polaron crystallizatio n 31 ' 32 into a 
metallic phase. Correlation-driven localization can be 
turned into correlation-driven superconductivity by small 
modification of certain tunable experimental parameters. 
Real systems where such and similar phase changes are 
manifest (cuprates, manganites, transition metal oxides, 
magnesium diborides, just to name a few) are complex 
in their structure as well as in the composition of their 
metal-ion-ligand electronic configurations. Their partic- 
ularity lies in a prevailing local physics, which has been 
generally overlooked in theoretical studies of the polaron 
problem. This local physics nevertheless dictates the out- 
come of macroscopic states such as (i) a superconducting 
phase of itinerant electrons with Cooper pairing, (ii) an 
insulating phase of localized bipolarons, and, possibly, 
(iii) a phase corresponding to a Bose metal sandwiched 
between the two, which has become of particular interest 
in recent years i 33 ' 34 ' 35 Since the physics of electron-lattice 
coupled systems involves local atomic or molecular dis- 
placements, it requires a description in terms of small 
polarons, i.e., electrons dragging with them a sizable lo- 
cal lattice deformation which can be of the order of a 
tenth or so of the lattice constant. We have for that rea- 
son to rely on the picture of a Holstein-type molecular 
crystal model as the basis of our study. 

As a paradigm for systems with resonant pairing in- 
duced by local lattice deformations, we can imagine as 
the simplest scenario a bipartite lattice like the one illus- 
trated in Fig. 1 and consisting of: 

(A) a sublattice of polaronic sites, where the charge 
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FIG. 1: Paradigm of a 2D resonant pairing system in terms of 
a bipartite lattice composed of polaronic cation-ligand complexes 
(rings) comprising four anions (large circles) in the center of inter- 
linked plaquettes of four metallic cation sites (small circles) housing 
the charge carriers in itinerant states. 



carriers couple locally and strongly to the molecular de- 
formations of the ligands, resulting in localized bipola- 
ronic bound pairs in the cation-ligand bonds; 

(B) a sublattice of non-polaronic sites on which the 
electrons (holes) move as itinerant charge carriers. 

When the charge transfer between the two sublattices 
is switched on, the electron (hole) pairing on polaronic 
sites, leading to localized bipolarons, will compete ener- 
getically with the itinerancy of uncorrelated charge car- 
riers on the non-polaronic sites in their immediate vicin- 
ity, when their respective energies are comparable. This 
induces pair correlations among the otherwise intrinsi- 
cally uncorrelated itinerant charge carriers in locally very 
confined regions of the lattice. On a macroscopic scale, 
this local physics (when treated on the basis of a heuris- 
tic Boson-Fermion Model scenario) ultimately results in 
a competitions between a correlation-driven insulating 
phase and a superconducting one i 35 ' 36 The transition 
from one phase into the other can be tuned by varying 
the strength of the electron-lattice coupling, the size of 
the adiabaticity ratio, the relative ionic potential differ- 
ence between the polaronic and the non-polaronic sites, 
and the concentration of the charge carriers. Specific in- 
terconnected local properties in such systems, which are 
determinant for the physics on a macroscopic scale, are: 

i) the generation of an intrinsic pseudogap feature in 
the local electronic density of state o 34 ' 35 which, depend- 
ing on the parameters of the system, foreshadows an in- 
sulating or a superconducting gap; 

ii) the strong softening and a concomitant increase of 
spectral weight of local vibrational modes on the pola- 



ronic sites^ serving as a basic ingredient for the bipo- 
larons to acquire itinerancy and assuring a component of 
bosonic charge carrier transport; 

(iii) the double-peak structure in the local lattice dis- 
placement pair distribution function due to a quantum 
coherent superposition 24 of bipolarons and pairs of intrin- 
sically uncorrelated electrons in their immediate vicinity, 
which can result in dynamically fluctuating stripe-type 
topological structures. 

The approach to the polaron problem followed here 
is similar to procedures followed in the correlation 
problem^ based on Anderson's Resonating Valence 
Bond idea^. There, an intrinsic underlying magnetic 
structure leads to the formation of local singlet hole 
pairs on plaquettes, dynamically exchanging with pairs 
of uncorrelated itinerant holes in their immediate vicin- 
ity. As concerns the connections with other electron- 
phonon models, the present paradigm of a fluctuating 
bipartite lattice structure of polaronic and non-polaronic 
effective sites is qualitatively different from the usually 
studied molecular crystal Holstein model, where all the 
sites are considered to be equivalent. Moreover, when 
the attention is focused on charge Kondo-like resonance 
featureS ) 39 i 40 ' 41 the local electron-phonon coupling is as- 
sumed to be of the form tuvoafoi^ + n i,l ~ ^)i a t + a «L 
i.e. it is at any site symmetric with respect to the single- 
site occupation. On the other hand, as we have stressed 
in the beginning, systems with double charge fluctuations 
show strong anisotropy of their ligand deformations. This 
motivates us to take an asymmetric electron-lattice cou- 
pling of the form Huoctfoi^ + n^i^af + <2j]. Although 
the modulus of the displacements of either unoccupied 
or doubly occupied sites relative to the singly occupied 
sites, i.e. (a^ + a,) = ±201(77^1+77^1), is the same for the 
two types of couplings (symmetric and asymmetric), the 
respective energies, related to the electron-lattice cou- 
pling, are not. For the usual symmetric electron-phonon 
coupling the energies of the unoccupied and doubly oc- 
cupied site are degenerate. It is this which results in the 
Kondo-type resonance in the charge channel ) 39 i 40 ' 41 pro- 
vided that the chemical potential of the system is such 
that this degeneracy, representing the two stable valence 
configurations, is guaranteed, which is the case for the 
exactly half-filled band limit. Away from this limit this 
is not the case and the resonance disappears A A fur- 
ther and essential difference between the generic Holstein 
scenario and the present paradigm for resonant pairing 
systems, is that in the valence fluctuation scenario the 
ligands of two adjacent cations have an anion in com- 
mon. This leads to strong intersite correlations between 
cations with different valence states and correspondingly 
different ligand environments compatible with this bridg- 
ing oxygen. These facts differentiate between sites where 
the bipolarons form and sites in the immediate vicinity 
where they do not. This is similar to what happens with 
Zhang-Rice singlets in cuprates ; 42 ' 43 where singlet pairs 
form on a Cu ion and its immediate ligand environment 
of four oxygens, but where nearest-neighboring Cu ions 
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cannot behave as such Zhang-Rice singlets at the same 
time. This "exclusion" feature^ introduces in the end 
strong local correlation effects in the band structure of 
the whole system, which is part of the salient features of 
the cuprates. Our proposition here is that the local pola- 
ronic features of double charge fluctuations should have 
a similar effect on the overall band structure and in par- 
ticular will give rise to a charge pseudogap of polaronic 
origin. 



III. THE MODEL 

The aim of this work is to examine the local physics 
of polaronic systems lending themselves to double charge 
fluctuations, coupled via an exchange term to itinerant 
charge carriers. They impose on the latter strong local 
dynamical diamagnetic fluctuations which can eventually 
lead to a global phase-locked superconducting state. We 
shall do this here, having in mind the cuprates and in 
particular their normal phase properties. These systems 
are more complicated than the paradigm for resonant 
pairing systems discussed in the previous Section, in the 
sense that all the Cu cations are intrinsically equivalent 
and there is no translational symmetry breaking to be ex- 
pected from the outset. Yet, we had early-on strong ex- 
perimental indications^ that the structural correlations 
and deformations in the cuprates show remnant bipartite 
lattice features, on a local level as well as on a finite time 
scale. These topologically different features, generally re- 
ferred to as " stripes" , break the translational, as well as 
the rotational, symmetry on a local level. 

The application of the valence fluctuation scenario, de- 
scribing possible resonant pairing in high-T c supercon- 
ducting cuprates, implies, for the reasons clarified below, 
the consideration of three formal valence configurations 
in the Cu02 layers, Cu 1 , Cu 11 and Cu 111 . Distinct lig- 
and deformations are associated with them, dumb-bell 
for Cu 1 and square planar for Cu 11 and Cu m , but dif- 
fering from each other by the value of the Cu— O dis- 
tances (see Fig. 2). The Cu cations in the CuC>2 layers 
of the insulating parent compound have a square oxygen 
ligand environment with Cu— O distances equal to 1.94 
A and are in formal valence states Cu 11 , according to 
their chemical composition. Moreover, the cuprates rep- 
resent a metastable structure characterized by a geomet- 
ric misfit of the various layers^ which leads to Cu— Cu 
distances in the Cu02 layers which can accommodate 
such [Cu 11 — O4] units (see Fig. 2a) only by pushing the 
oxygens alternately slightly above and below these lay- 
ers. This structural effect is usually called buckling. 

Upon hole doping one introduces single holes in such 
[Cu 11 — O4] units, which make them turn into [Cu m — O4] 
units, with shortened Cu— O distances equal to 1.84 A. 
Referring to our small cluster, such a doped effective 
cation-ligand unit is represented by the central effective 
site 5 in Fig. 2b. The reduction to 1.84 A of the bond- 
lengths in the [Cu 111 — O4] units is accompanied by an 



unbuckling of the four oxygen cations belonging to these 
units, which arc pulled into the Cu02 plane, but keep 
unaltered their distances from the Cu 11 ions on the pla- 
quette sites, thus inducing no long-range stress field. In 
this way, the oxygen cations surrounding each of the Cu 11 
ions (located at sites 1-4 in Fig. 2b) are still all at the 
same distance 1.94 A from it, but one of them [the bridg- 
ing oxygen between CU5 and Cu; (i=l — 4)] now lies in 
the Cu02 plane and not above or below it. 

We conjecture that the [Cu 111 — O4] units with their 
doping-induced ligand deformations act as polaronic 
traps, binding two holes in form of a localized bipolaron. 
If the energy of this state is almost degenerate with that 
of the hole states at the Fermi level, then this bipolaron 
can split into two holes going on the surrounding matrix, 
and, as a consequence, [Cu 111 — O4] turns into [Cu 1 — O4]. 
As pointed out in the previous Section, Cu 1 is typically 
in a linear O— Cu— O dumb-bell configuration. If we sup- 
pose that this happens also in the Cu02 layers, we are 
faced with two degenerate O— Cu— O dumbbell bonds 
along the x and y directions, with a shortened bond- 
length of 1.84 A (see Fig. 2c) which is the same as for the 
[Cu m — O4] units. This would cause an anisotropic un- 
buckling of the corresponding bridging oxygens in cither 
one or the other direction, but leaving the Cu 11 — O4 units 
on the plaquette sites unaltered, as far as their Cu— O 
bond-lengths are concerned. It is also feasible that, given 
this directional degeneracy of the bonds of the Cu 1 va- 
lence states, an isotropic ligand deformation with the 
same bond-length of 1.84 A materializes, due to d x 2_ y 2 
hybridization in the CuC>2 planed and which then would 
be very similar to the configuration illustrated for the 
Cu 111 — O4 unit in Fig. 2b. It would lead to an isotropic 
unbuckling, again with no changes in the bond-lengths of 
the cation-ligand units on the plaquettes. 

The microscopic foundation of the various formal va- 
lence configurations, conjectured here to play a role in 
the cuprates, are left to be investigated in the future. 
For the present study, where we shall treat the lig- 
ands as an effective deformable scalar quantity, we do 
not deal with such detailed questions concerning the 
anisotropy of the deformations on the central sites act- 
ing as resonating bipolaronic traps. The valence fluc- 
tuations involving the central [Cu 111 — 04)5 unit lead to 
the release of two holes from the latter onto the neigh- 
boring [Cu 11 — O4]; units on the plaquette sites i=l — 4, 
and vice- versa ([Cu 111 — 04)5 ^-[Cu 1 — O4] 5 + two holes). 
Once on the plaquette units, these holes engage in itin- 
erant states which prevent them from getting localized 
in form of self-trapped bipolarons on those [Cu 11 — O4]; 
units. Such a mechanism describes tunnelling in and 
out of hole pairs from a given central deformable ligand 
[Cu— 04)5 unit (our effective site 5) into the surrounding 
[Cu 11 — O^j units which remain essentially undeformed. 
For this reason we neglect here any charge-ligand defor- 
mation coupling on those plaquette sites. 

We remind the reader that formal valence in such sys- 
tems has to be interpreted in terms of the charges in- 



volved in the covalent bonds between the cations and 
their ligand environments and not simply by the charges 
on the cation. In the cuprates the covalent character in 
the C11O2 planes is very pronounced. This might be one 
reason why spectroscopically neither Cu 111 nor Cu 1 have 
been detected sofar. Furthermore, given the fact that 
the three formal valence states Cu 1 , Cu 11 and Cu 111 have 
different ligand environments, the coupling of those de- 
formations to the charge will in general not be symmetric 
with respect to singly occupied units [Cu 11 — 4 ]. In the 
planar Cu02 layer of the cuprates all the Cu ions are 
equivalent and their charge coupling to the ligand defor- 
mations are spatially homogeneous. But this coupling 
does depend on the formal Cu valence state which is ac- 
tive at any given moment. If we consider the cuprates 
as being made out of an ensemble of overlapping clus- 
ters, as given in Fig. 2, this immediately leads to local 
dynamical correlations between neighboring ligand defor- 
mations. On a short range we can therefore expect that 
this will naturally lead to intercalated topological local 
structures. The anisotropy of the O— Cui— O bonds could 
be relevant for the anisotropic local symmetry breaking. 

The construction of a macroscopic state on an infi- 
nite lattice, starting from a spatially translationally in- 
variant Hamiltonian but working in a representation of 
such finite cluster states, which contain the essential local 
physics, can in principle be achieved. One way would be 
using renormalization group techniques such as the pla- 
quette contractor method, based on small clusters of the 
form illustrated in Fig. 2. This method was developed 
by Morningstar and Weinsteini^ and applied to the RVB 
scenario by Altman and Auerbach 37 . For the present 
scenario of valence fluctuation driven pairing, including 
the various ligand deformations correlated to the formal 
valence state of the clusters, presents a formidable nu- 
merical task. Questions like that will be addressed in 
some future work. 

Here we shall content ourselves to study the dynam- 
ics of such basic clusters making up the Cu02 layers, 
which permits us to draw some preliminary conclusions 
on such macroscopic phases. Considering an ensemble of 
such spatially and temporarily uncorrelated fluctuating 
clusters, without any long-range spatial phase coherence 
between them, allows us to describe certain features of 
the pseudogap phase in the cuprates and its correspond- 
ing onset temperature T* . We can assess the frequency 
of the resonant pair tunneling processes, driven by the lo- 
cal ligand deformation fluctuations (which play the role 
of local dynamical lattice instabilities in a macroscopic 
system) and which can be approximately related to the 
mass of the diamagnetic pair fluctuations and thus to the 
value of T c for the onset of a Bose- Einstein condensation- 
driven superconductivity. 

In order to capture the salient features of the local 
physics of valence fluctuations in the cuprates, as exposed 
above, we illustrate in Figs. 2a-c. the various valence con- 
figurations of such clusters, where the central Cu cation 
is alternatively in a Cu 1 , Cu 11 and Cu 111 formal valence 
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FIG. 2: (Color online) Basic cluster in double valence fluctuation- 
driven pairing in the CUO2 planes consisting a central Cu cation 
oxygen ligand environment in three different formal valence states, 
Cu 11 , Cu 111 and Cu 1 , illustrated in Figs. 2a, 2b and 2c, respectively. 
Such central cation-ligand units are linked via a common bridging 
oxygen to neighboring Cu n -04 units. The small red circles denote 
the Cu cations and the large circles the oxygens. Dark blue means 
that they are pushed above the CUO2 planes, light blue that they 
are below and white circles that they are in-plane. 



state. We address this local problem here in a some- 
what simplified version, neglecting the detailed structure 
of the oxygen ligand environments surrounding the Cu 
cations. For that purpose we assume the [Cu— O4J5 unit 
to be composed of (i) a central cation surrounded by 
an isotropically deformable ligand environment of four 
anions on an effective site 5, and (ii) four neighboring 
[Cu— O4]; units (i=l-4) on the plaquette, representing 
non-polaronic sites for which we consider the charge- 
ligand deformation coupling to be inactive. This is be- 
cause we consider the charge carriers on the plaquette 
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sites to move as itinerant entities and those ligands have 
no time to respond to the charge-deformation coupling. 
They remain in their buckled configurations, correspond- 
ing to the [Cu 11 — 04]; units. We shall denote the hopping 
rate between the plaquette sites by t and that between 
the latter and the central polaronic deformable unit by 
t* . The dimensionless coupling between the charge den- 
sity and the local ligand deformations is denoted by a 
and the bare vibrational frequency of the bond-length 
fluctuations of that latter by uj . We assume that the 
ionic levels of the plaquette sites and of the polaronic 
sites have a difference in energy, given by A. The reason 
for that is the following. Our picture is that the local 
clusters, as illustrated in Fig. 2, represent the essential 
dynamical units of the cuprate Cu02 planes. In this pic- 
ture, the states of holes on the plaquette sites play the 
role of the states of itinerant holes in an infinite system 
close to the Fermi level, moving in a Hubbard correlated 
fashion. The holes we tackle are those induced by dop- 
ing the parent compound of the correlated problem but 
which, on an individual cluster, we shall simply treat as 
itinerant uncorrelated holes on the plaquette. Includ- 
ing the Hubbard U in this problem would require one to 
consider at the same time the electronic correlations on 
the polaronic sites, leading to Zhang-Rice singlets, and 
those on the plaquette sites, reflecting their pairing in 
the overall density of states near the Fermi level. This is, 
however, beyond the present description. 

Given these considerations, the Hamiltonian describing 
such cluster is 



H 



* E 

i=tj=l..A,<j 



c ia C 3,<r + h.C. 



E [ C V c 5,v + h -c] + A E C L C 5, 



a) a + 



- huj a ^ c l,a c 5,a [a + a f ] , (1) 



where c^J denote the annihilation (creation) operator for 

a hole with spin a on site i, and a 5 the phonon annihila- 
tion (creation) operator associated with the deformable 
cation-ligand complex at site 5. This local Hamilto- 
nian describes a competition between localized bipola- 
ronic hole pairs on the central cation-ligand complex and 
itinerant holes on the plaquette sites, when the ener- 
gies of the two configurations are comparable, i.e. for 
2A — AfvjjQa 2 ~ — At. We illustrate in Fig. 3 the varia- 
tion with a of the low energy spectrum of such a system 
for a set of parameters (t = 0.2, t* — 0.15, A = 0.5 
and ojq — 0.1, in bare units) corresponding to a typical 
region in the parameter space [u>o/t, hujoa 2 /t] where reso- 
nant pairing is well pronounced and rather robust against 
thermal fluctuations. Our approach is a straightforward 
exact diagonalization of H leading to eigenenergies E n 
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FIG. 3: (Color online) The low energy two-hole spectrum of H 
as a function of the electron-lattice coupling constant a, showing 
the crossover between delocalized itinerant states on the plaquette 
(for a < ct c = 2.12) and localized bipolaron states on the central 
polaronic site (for a > a c ). 



associated with eigenstates of the form 



n=EE4>)®h 



(2) 



k u=0 



where designates the configurations of two holes with 
opposite spin distributed over the five cluster sites, \v) 
denotes the ^-th excited state of the undisplaced har- 
monic oscillator state, and A^l is the corresponding 
weight in the eigenstate \n). Taking into account a trun- 
cated Hilbert space for the phonon states \v), we limit 
ourselves to 7V p h = 70. For the set of parameters and 
temperatures we are concerned with, this is sufficient. 
We notice from Fig. 3 that the low energy eigenvalues 
E n ~ — At + nhujQ are relatively unaffected as a is in- 
creased, except when they get close to crossing the en- 
ergy of the localized bipolaron level at some a c . As a is 
increased above this value, the energy eigenvalues follow 
the downward trend of the localized bipolaron energy, 
tending to E n ~ 2A — Ahw^a 2 + nhu>o. 

In Fig. 4 we trace out the occupation probabili- 
ties in the two-particle ground state \GS) for the 
configurations with (i) the two holes on the pola- 
ronic site, (GS\c' 5 i-cSj-fCg iCs^ilGS), (ii) the two holes 
on the plaquette sites, (GS'|cJ =0 jC g= o,T<4 =0 i c q=°>l\GS)i 
where c^ =Ocr |0) = | C ! CT |0) signifies the lowest- 
energy Bloch state of the plaquette, and (iii) the 
two holes distributed between the two types of sites, 



E CT =-r,i( GS 'l c 5 ) <7 C 5,<rcJ = o_ .Cg=o,- CT |GS'). Fig. 4 shows a 
rapid change-over when we sweep a through a narrow 
regime around a c = 2.12 for u>o/t = 0.5. It is in this nar- 
row regime that we can expect resonant pairing, which 
means a tunneling of bipolaronic hole pairs from the cen- 
tral site to the plaquette sites and vice- versa, in this way 
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FIG. 4: (Color online) The occupation probabilities for the various 
possible hole configurations as a function of a and for uio/t = 0.5. 
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FIG. 5: (Color online) The efficiency factor FfLu as a function of 
bjQO? jt for a hole tunneling in and out of a polaronic site, favoring 
double charge fluctuations, for several adiabaticity ratios LUo/t- For 
clarity of presentation, the curves obtained for uJo/t = 0.5,1,1.5 
have been rigidly shifted upwards with respect to the one referring 
to uo/t = 0.375 by amounts equal to 0.04, 0.08, 0.12, respectively. 
The numerical values on the j/-axis thus refer specifically to the 
case u>o/t = 0.375. 



inducing a hole pairing on the latter. 

The efficiency for double charge fluctuations is con- 
trolled by the ease with which holes tunnel in and out of 
sites. At zero temperature it is given by the ground-state 
static correlator 

^cxch = (GS\ C q=0,<T C 5, cr\GS) (3) 

where \GS) denotes the ground state for the two-hole 
system. The variation of F^ ch with a is illustrated in 
Fig. 5. We notice that for a close to zero, the efficiency 
for such tunneling is governed by the residual tunneling 
rate determined by the difference in energy between the 
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FIG. 6: (Color online) The efficiency factor Fl^-a. as a function 
of cooa 2 /t for converting, via resonant tunneling, a localized bipo- 
laron on the central polaronic cation-ligand complex into a pair of 
correlated holes on the ligand environment, for several adiabaticity 
ratios LJo/t. 



central site and the lower-lying Bloch states on the pla- 
quette, equal to A — 2t. As one approaches a certain 
a c from below, we notice a strong increase in this tun- 
neling efficiency. However, when going above a c , this 
efficiency drops rapidly to zero, meaning that the holes 
are now quasi-localized on the polaronic site in form of 
bipolarons. 

In order to have a double charge fluctuation-induced 
superconducting phase, driven by dynamical local lat- 
tice instabilities, the above criterion for a high efficiency 
rate of hole tunneling in and out of polaronic sites is not 
sufficient. What is relevant is the efficiency for coherent 
double charge fluctuations between the polaronic site and 
the plaquette. Such coherent, or resonant, pair tunneling 
is given by the efficiency factor 

^exch = (GS\4 iT 4 aC , = o,ic,=o,T \GS) [F c c ' ch ] 2 , (4) 

where we have subtracted out the rate of exchange of 
the bipolaron with two holes on the plaquette via inco- 
herent single-hole processes, given by the second term in 
this expression. The efficiency factor for resonant pair 
tunneling (see Fig. 6) is sharply peaked near a c , its max- 
imal intensity depending on the adiabaticity ratio u>o/t. 
The coherent pair tunneling represents only a small frac- 
tion of the total tunneling rate, determined largely by its 
incoherent part. Nevertheless, the size of this coherent 
component of the two-particle spectral function is suf- 
ficient to ensure a high enough condensate fraction in 
order for a superconducting phase-correlated state to be 
stabilized. It provides a strong support in favor for such 
a mechanism to play a role in real materials which have 
tendencies to local dynamical lattice instabilities driving 
double valence fluctuations. 

In order to visualize resonant pair tunneling in and out 
of the polaronic cation-ligand complex, we examine the 
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FIG. 7: Evolution in time r of the resonant tunneling of a hole 
pair in and out of a polaronic site, composed of a cation together 
with its surrounding ligand environment, for LOg/t = 0.5. 




FIG. 8: Evolution in time r of the bond-length of the ligand envi- 
ronment which accompanies the resonant tunneling of a hole pair 
in and out of the cation-ligand complex, for uJo/t = 0.5. 



evolution with time r of double charge fluctuations (see 
Fig. 7). It is described by the correlation function 



(GS\[n k =o,t(r)n k =o,l( T ) - n 5;T (r)n 5a (T)] 
K = o, T (OH =(U (0) - n 5;T (0)n 5li (0)]G5). 



(5) 



Near a c = 2.12 (for Lo /t = 0.5) the charge disproportion- 
ation follows a smooth oscillation in time with a period 
ujqt ~ 70 and an amplitude varying between +0.5 and 
—0.5. When this correlation function is positive and at 
its maximum, such as for uj t = 0,70,..., it indicates 
that both holes are on the central polaronic site. If it is 
negative and maximal, such as for loqt = 35,105,..., it 
means that both holes are on the plaquette. The tun- 
neling frequency of the coherent part of this motion of 
the hole pairs is u>q = [2ir / loot]luo — [27r/70]wo — 0.09o>0j 
as we see from Fig. 7. It represents the frequency for 
a correlated mode, locking together the dynamics of 
the charge and deformations of the ligand environment, 
or, in our simplified description, of some bond-length 



X = y / h/2Muj [a^ + a] . The fluctuations of the latter 
are described by the correlator 



C fl b u ° c n t d (r) = (GS|[ot(r) + a(r)][ O t(0) + a(0)}\GS) (6) 

and its variation in time r is illustrated in Fig. 8. In 
line with the fluctuations of the charge distribution, 
this bond-length fluctuation shows for a ~ a c the 
same smooth oscillating behavior with frequency Wq — 
0.09, centered around some mean deformation (GS'lfcr + 
a] ICS*}. The important feature of such resonant pair tun- 
neling is that the fluctuations, which accompany the sys- 
tematic motion of the charge distribution as well as of 
the bond-length, 

(i) arc on a time scale, typically of the order of l/^o, 
much shorter than the one of the coherent tunneling mo- 
tion ~ 1/wq; 

(ii) have an amplitude which is small compared to the 
amplitude of the systematic coherent oscillating behavior 
on the time scale of l /ojq . 

As soon as we move away from this resonance regime, 
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by tuning a away from a c for a given adiabaticity ra- 
tio, the bond-length fluctuations very quickly lose their 
smooth coherent oscillating behavior. They begin to 
show phase slips, which, as a result, decorrelate the hole 
pair fluctuations from those of the bond-length ones. As 
a result, the hole pairs get confined on either the central 
polaronic site for a > 2.14 or else on the plaquette sites 
for a < 2.10. In our small system this confinement is in- 
ferred from the absence of any sign of periodic motion of 
the charge and bond-length fluctuations between those 
two different kinds of sites. This is true for the whole 
regime of ljqt between zero up to at least 20000. For all 
intents and purposes, one can therefore regard such a be- 
havior either as a localization for a > a c of tightly bound 
hole pairs on the polaronic site, or as a confinement for 
a < a c of two uncorrelated holes on the plaquette. In- 
corporating this feature into an infinite lattice structure, 
such as the one depicted in Fig. 1, we expect that for 
a > a c we have a true charge disproportionation, i.e., 
the system breaks up into a mixture of cations with lo- 
calized bipolaronic hole pairs on the cation-ligand sites 
and leaving the remaining cations without any hole. For 
the opposite case, a < a c , we expect a metallic behavior 
with locally correlated pairing of holes induced by a dy- 
namical charge exchange coupling with the neighboring 
polaronic sites. Depending on the overall particle den- 
sity, this can give rise to: (i) charge-ordered states with 
translational symmetry breaking, passing from a cova- 
lent to ionic bonding^ (ii) correlation-driven insulating 



i 1 1 1 r 




IV. SIGNATURES OF RESONATING 
BIPOLARONIC DOUBLE CHARGE 
FLUCTUATIONS 

The phase-correlated dynamics of the charge and the 
ligand deformation results in several characteristic fea- 
tures in the electronic and lattice properties which are in 
principle tractable experimentally. Since these features 
arise from purely local correlations, they are adequately 
described by the small cluster model system studied here. 
These features involve: 

(i) The strong softening of a local phonon mode with 
frequency u>q (see top panel of Fig. 9) which describes 
correlated fluctuations of the charge and the local lattice 
deformations of the cation-ligand complex. This mode is 
described by the spectral function for bond-length fluc- 
tuations, given by 

SphM = --lmDph(w) (7) 

with 

£> ph (o;) = ^°((X;X)) w . (8) 



FIG. 9: (Color online) Lattice dynamical signatures related to 
resonant bipolaronic double charge fluctuations for several adia- 
baticity ratios u)$/t and as a function of the electron-lattice cou- 
pling ujqoi 2 /t. Top panel shows the softening of the intrinsic bond- 
length vibrational mode as one approaches the critical coupling 
for resonance tunneling at a c ~ \J[t + A/2)/uo, depending on 
ujo/t. The maximal softening, occurring for ujo/t = 0.375, 0.5, 1, 1.5 
at a c = 2.45,2.12,1.5,1.22, respectively, is given by U!q (a c )/t = 
w*/t = 0.00336, 0.00907, 0.04205, 0.08103. The central panel shows 
the associated evolution of the spectral weight £? p i 1 (w = u)q) of this 
renormalized mode. The bottom panel shows the spectral weight 
Bpi^u! = 0) of the zero-frequency mode, related to a static overall 
shift of the bond-length. 



At finite temperature i? p h(w) is given by 

^phM = ~^e-^\(m\(a* + a)\n)\ 2 5(u + E n -E m ) 
n.rn 

(9) 

with Z = J2 n e ■ A s this local mode softens upon ap- 
proaching a c , its spectral weight Bp^us = Wr?) strongly 
increases, reaching a high maximum at uiq(ol c ) = ujq (see 
middle panel of Fig. 9). This is indicative of the large 
amplitude fluctuations of the bond-length which are nec- 
essary to drive the double charge fluctuations between 
the polaronic cation-ligand complex and the plaquette 
sites in a coherent fashion. 

(ii) The spectral weight of the zero frequency bond- 
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length correlation function Bp^ui — 0) is practically zero 
for a < a c , sharply increases in a narrow region around 
a c , and then continues to grow linearly upon further in- 
creasing a (see bottom panel of Fig. 9). This is indica- 
tive of a dynamical local lattice instability which should 
manifest itself as a central peak structure - a quasi-elastic 
peak - in the neutron scattering cross-section. Such 
features have been observed previously in systems with 
local dynamical lattice instabilities near ferro-electric 
phase transitions^ martensitic phase transition in A15 
compounds^ and SrTi0 3 i£2i^ 

(iii) The emergence of a double-peak structure in the 
static Pair Distribution Function (PDF) for the bond- 
length, 

n 

= |E e ^"EE4^>i^-w> 

n k v % v' 

= E 7nP0 (10) 
n 

which can be measured by neutron spectroscopy or ex- 
tended X-ray absorption fine structure (EXAFS). This 
function has been investigated in experimental studies 
concerning the local dynamical lattice topology in the 
cuprate superconductors and in the manganites, for the 
Cu— O and for the Mn— O bond stretch mode, respec- 
tively (for a review see the contributions by Sinai et al. 
and by Egami in Ref. 25). 

Our results, shown in Fig. 10, indicate how the PDF 
evolves near the critical value of a c , where resonant pair- 
ing is present and how it is affected by the temperature. 
The left column of this figure shows its evolution as a 
crosses a c , for several temperatures. At T = there is a 
clear signature for a double-peak structure at a c , quickly 
disappearing as we go away from this resonance condi- 
tion. At very low temperatures the double-peak struc- 
tured PDF can be taken as an indication for resonant 
pairing. At finite temperature, however, one cannot ac- 
credit the double peak to a resonant tunneling of double 
charge fluctuations. This becomes clear from the second 
and third column of this figure, showing the contribu- 
tions 7o(AT) and 71 (X) to the total PDF, coming from the 
ground state and the first excited state, respectively. At 
finite temperature, a double-peak structure always arises 
away from a c from thermal excitations which contribute 
to the PDF in form of phase-uncorrelated single-peak 
contributions. 

(iv) The onset of strong local pair correlations between 
the charge carriers below a certain temperature T*, re- 
sulting in a pseudogap in the local density of states of 
our cluster. It was in fact on the basis of the present sce- 
nario that a pseudogap phase was initially predicted for 
the cuprates above the superconducting transitioni 55 i 56 
Its lattice-driven origin became apparent in the isotope 
effect experimentally observed 5 ^ on the temperature T* 



at which the pseudogap opens up. 58 Within the same sce- 
nario, the onset of the pseudogap feature was found to be 
related to a change from single-particle transport at tem- 
perature above the onset of the pseudogap (T > T*) to 
one which is controlled primarily by bound charge carrier 
pairs below T* , upon approaching T c from above. These 
bound pairs then acquire free-particle-like dispersion, 
which results in a transient Meissner effect 5 ^ (experimen- 
tally verified^) and remnant Bogoliubov shadow modes 
in the pseudogap phased These theoretical results were 
obtained by assuming an effective local static exchange 
mechanism between bipolaronic bound hole pairs and 
itinerant uncorrelated pairs of holes. The present work 
shows to what extent such a heuristic scenario, investi- 
gated within the Boson-Fermion Model, holds true and 
to what extent a time dependent relaxation of the cation- 
ligand environment, trapping momentarily the charge 
carriers in form of bipolarons, permits one to use such 
an effective double-charge exchange mechanism. 

In Fig. 11 we illustrate separately the contributions to 
the DOS coming from the holes on the plaquette and 
those coming from the central polaronic site. The first 
ones involve the region near the center of the DOS and 
show a significant amplitude as well as asymmetry. That 
latter is a consequence of the asymmetric charge-lattice 
deformation coupling which we assumed in this work (see 
Section 2). The holes on the polaronic site are respon- 
sible for the peaks on either side of the center of the 
DOS at large frequencies. We notice that as the tem- 
perature decreases (right hand panels of this figure) the 
pseudogap-like structure turns into a well defined gap. It 
is this feature which, on the basis of an effective Boson- 
Fermion Model treatment ) 34 ' 35 finally leads to a true su- 
perconducting gap below T c . 

(v) Resonant pair tunneling in and out of the polaronic 
cation-ligand complex into the surrounding sites on the 
plaquette incites pair correlations in the latter. They 
give rise to the residual pairing in a many-body state, 
where local diamagnetic fluctuations are embedded in an 
underlying Fermi sea and eventually result in a macro- 
scopic superconducting phase. Indications for such local 
diamagnetism in the pseudogap phase of the cuprates 
have been observed^, We illustrate in Fig. 12 the varia- 
tion with a (near a c ) of the average density n p 

n p = (n q=QA n q=QA ) - (n q=0yl )(n q=od ) (11) 

of diamagnetic paired holes on the plaquette and notice 
that is it quite sizeable - about 20% of the total aver- 
age density of holes on the plaquette at a c . Slightly 
away from this resonant regime, the holes are either 
pair-uncorrelated (for a < a c ) or not existent on the 
plaquette (for a > a c ), being confined to the polaronic 
cation-ligand complex as localized bipolarons. The vari- 
ation with the temperature of n p is illustrated in Fig. 13. 
The contributions to n p come primarily from two sources. 
The first and dominant one comes from the first excited 
state of the system, lying an energy ujq above the ground 
state, and reflects the resonant nature of this diamagnetic 
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FIG. 10: (Color online) Left panels: pair distribution function, 
g(X), measuring the instantaneous distribution of the cation-ligand 
bond-length at the central site, for three values of a near a c and 
three different temperatures, for uio/t = 0.5. Central panels: pair 
distribution function at T = 0.05, together with the partial contri- 
butions 70 (X) and 71 (X) coming from the ground state and from 
the first excited state, respectively. Right panels: same as in the 
central panels, for T = 0.1. 
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pairing. The second contribution comes from the ground 
state itself, which indicates a static diamagnetic polariza- 
tion of the holes induced by this resonant pairing effect. 
As a result, as the temperature increases, n p initially in- 
creases, then, it flattens off at a certain temperature T*, 
and eventually decreases to zero due to thermal fluctu- 
ations. We thus take T* as a measure of the tempera- 
ture for the onset of resonant pairing and the opening of 
the pseudogap of the electronic DOS. The diamagnetism 
is strongest in the adiabatic regime where the contribu- 
tion coming from the resonant pair tunneling is maxi- 
mal, showing a peak of n p at low temperatures. These 
resonant features are clearly reflected in the dynamical 
pairing susceptibility (not presented here) which shows 
a strongly peaked behavior at a frequency uj^ above the 
ground state energy, indicative of the resonant nature of 
such diamagnetic fluctuations. 



V. SUMMARY 

The main objective of the present work was to exam- 
ine to what extent double charge fluctuations, driven by 
local dynamical lattice instabilities, can result in reso- 
nant pair tunneling in and out of polaronic cation-ligand 
complexes, such that spatial phase coherence of these 
fluctuations on a long-range scale can ultimately be en- 
visaged. This local physics with coherent resonant pair 
tunneling is a prerequisite for ensuring the existence of 
itinerant bosonic quasi-particles describing diamagnetic 
local pair fluctuations, embedded in a Fermi sea. It re- 



FIG. 11: (Color online) Single-electron DOS decomposed into 
the contributions coming from the plaquette sites (the spectral 
lines near the center of the DOS, in red) and from the central 
polaronic site (the spectral lines at the wings at high frequen- 
cies, in green), (a) Evolution of the DOS at low temperature 
T = 0.05 (left column) as a crosses the critical value for reso- 
nant pairing at a = a c = 2.12. (b) Evolution of the DOS with 
increasing temperature at a c = 2.12 for uio/i = 0.5 (right col- 
umn), which shows the filling in of the pseudogap with increasing 
T. The energy u> is in bare units and the pseudogap evolves around 
A - 2hw a 2 c ~ -2t = -0.4. 



suits in the opening of a pseudogap in the single-particle 
DOS at T*, which above T c simulates that of a semi- 
conducting phase at finite T. As T is decreased, this 
pseudogap goes over smoothly into the superconducting 
gap below T c where part of it is still controlled by the 
gap arising from pair correlations in the normal state. 35 
In the strong coupling anti-adiabatic regime, resonant 
pair tunneling would reduce to an exclusively bosonic 
system, composed of localized bipolarons forming static 
charge disproportionation insulators^ Superconductiv- 
ity in systems with resonant pairing is qualitatively dif- 
ferent from that of a BCS state in the sense that its T c is 
controlled by phase rather than amplitude fluctuations, 
as is the case for the cuprate high-T c superconductors 
in the underdoped regime. This requires, on a dynami- 
cal basis, a correlated mixture of diamagnetic pairs and 
quasi-free charge carriers, as described by the Boson- 
Fermion Model in its simplest version, where Hubbard 
correlation effects, which are certainly important in the 
cuprates, are neglected. 

Using the values of and n p deduced from the nu- 
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2 2.5 

co a 2 /t 



FIG. 12: (Color online) Density n p of correlated hole pairs on 
the plaquette as a function of ujQa 2 /t, showing a strong enhance- 
ment close to a c ~ ^/ ((t + A/2)/uo due to resonant pair tun- 
neling. At a c , rip is typically of the order of 0.1, which has to 
be compared with the average total density of holes on the plaque- 
tte, 5Z CT ('-''S1 n 9=0,T' 1 q=0,J.}! which from Fig. 4 follows to be roughly 
equal to 0.5. 




FIG. 13: (Color online) Variation with the temperature T of the 
density n p of correlated hole pairs on the plaquette for various adi- 
abaticity ratios u>o/t and the corresponding o c 's. With increasing 
T, n p rises and then flattens off at some T*, which indicates the 
onset temperature for pair correlations on the plaquettes. 



merical results illustrated in Figs. 9 and 12, and referring 
to the adiabaticity ratios u>o/t = 0.375,0.5,1.0,1.5, we 
estimate for the superconducting transition temperature 
the values T c ~ (£/a)[4, 10, 20, 27] K. Here we have used 
the approximate expression T c ~ h 2 (^/a)(n p /m p a 2 ) for 
phase fluctuation-driven superconductivity, where £ de- 
notes the coherence length, a the lattice constant, n p and 
m p the number and the mass of the diamagnetic pairs. 
Trip is related to the resonant pair tunneling frequency 



by h 2 /2m p 



OJna 



2 . Within our scenario of a lattice 



instability-driven superconductivity, already such a crude 
estimate gives values for T c which are quite reasonable. 
They have been obtained by using a single input parame- 
ter, luq = 50 meV, taken from experiments concerning the 
Cu— O bond stretch mode. It is this mode which shows 
an anomalous softening for large g-vectors close to the 
zone boundary and is thought to play an active role in 
such possible dynamical local lattice instabilities in the 
cuprate superconductors! 63 ' 64 The value of the electron- 
phonon coupling a in this estimate for T c is taken at 
the resonance, i.e., a — a c = [2.45,2.12, 1.5, 1.22] for the 
various values of the adiabaticity ratios considered here. 
The nearest-neighbor and next-nearest-neighbor hopping 
integrals are taken to be of the same order of magnitude, 
i.e., t/t* = 0.2/0.15. 

We can also estimate the onset temperature T* of pair- 
ing correlations in the itinerant subsector of the charge 
carriers. From the behavior of the pairing amplitude n p , 
given in Fig. 13 as a function of the temperature T, T* 
should be approximately identified by the crossing point 
of the tangent to the rising portion of n p at low T and 
the horizontal line passing through the maximum of n p . 
This T* is sensitive to pressure which increases the lo- 
cal phonon frequency. Using the thermodynamic relation 
between pressure-induced changes ST* in T* and Su>q in 
the local phonon frequency Wq> 6T*/T* = (3(5ujq)/u)q 
(/? denoting a numerical factor of order unity, related to 
the Gruneisen parameter), this experimentally verified 
relation^ 5 - for the cuprates with (3 ~ — 3 holds for our set 
of values for u>q and T* with a f3 of order of unity. 

Our concept of resonant pairing-induced diamagnetic 
fluctuations bears some resemblance to what has been 
termed Cooper pairing due to over-screening. Over- 
screening the Coulomb interaction between the elec- 
trons implies a negative electronic contribution (besides 
the positive ionic one) to the static dielectric constant, 
e c \(oj = 0,q) < for certain q = k — k' vectors in 
the Brillouin zone and provided that it is smaller than 
— 1. This over-screening does not result in static global 
lattice instabilities^ 6 - It can thus provide a mechanism 
for superconducting pairing among holes with wave vec- 
tors [k, —k'] for much higher electron-phonon couplings 
than initially thought. This idea has been brought up 
again in connections with the cuprates superconductors, 
where the softening of the Cu— O bond stretch mode in 
the Brillouin zone along the directions of the anti-nodal 
wave-vectors, such as [0, ±0.5]7r/a and [±0.5, 0]n/a, was 
interpreted as a signature for over-screening, result- 
ing from local phase-correlated charge and bond-length 
fluctuations^ As we could confirm in the present work, 
such correlated pair fluctuations do indeed exist and are 
quite robust. However, based on our study, we do not ex- 
pect that they give rise to a BCS-type superconducting 
phase which is controlled by a pairing amplitude fluctu- 
ations. Such systems should be close to an insulator-to- 
superconductor transition ; 35 ' 36 the insulator being rep- 
resented by phase-uncorrelated pair fluctuations and the 
superconductor by spatially phase-correlated ones with 
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a T c controlled by phase and not amplitude fluctua- 
tions. As far as superconductivity is concerned, our 
approach to this proble m 34 ' 35 ! 36 ' 56 ! 58 differs from any 
standard Eliashberg formulation! 6 - 3 . In the Boson- Fermion 
Model scenario superconductivity comes about via a cen- 
ter of mass motion of dynamically fluctuating pairs. This 
is clearly demonstrated by (i) the early Uemura muon 
spin resonance experiments^ showing a T c varying with 
the concentration of superfluid carriers n p , and (ii) a 
T c ~ 2(£ / a)[n p ujQ] which decreases with increasing soft- 
ening of the bond-stretch mode frequency 6 -^ ujq — u)q, 
where ujq determines the inverse mass of these itinerant 
diamagnetic pair fluctuations. This experimental results- 
relates the doping dependence of T c to the local lattice 
properties, characterized by a softened local mode result- 
ing from dynamical local lattice instability for particular 
dopings in the superconducting regime. 

This tempts us to make the following conjecture for our 
resonating pair scenario, applied to the superconducting 
cuprates: Each doping level in these superconducting 
compounds correspond to the thermodynamically most 
unstable single-phase composition in the synthesis. Once 
stabilized in their solid phase at low temperature, the ki- 
netic stability of those metastable phases is determined 
by our conditions for resonant pairing. This fixes, for 
any given bare local mode with frequency ujq (relatively 
insensitive to doping), a critical charge- ligand deforma- 
tion coupling a c , which will depend on doping and con- 
sequently on the degree of softening of this mode. 

Resonant pair tunneling poses rather stringent and ex- 
perimentally verifiable conditions on the local physics, as 
discussed here. They are primarily related to features 
linked to dynamical lattice instability driven by charge 
fluctuations, such as: (i) the strong softening of a local 
phonon mode describing the correlated charge-ligand de- 
formation fluctuations (with a frequency ujq well below 
the bare local phonon frequency luq), (ii) a quasi-elastic 
peak in the neutron scattering cross-section indicative of 
the dynamical nature of a local lattice instability, (iii) 



a double-peaked pair distribution function, accessible to 
neutron scattering spectroscopy and EXAFS, (iv) locally 
fluctuating diamagnetism, and (v) a characteristic struc- 
ture of the electron spectral function, which shows po- 
laronic peaks at large frequencies on either side of the 
centre of the pseudogap and a strongly peaked DOS near 
this center, coming from the local diamagnetic pair cor- 
relations. 

The natural next step of the present analysis will have 
to incorporate these local dynamical features in a large 
scale system. Here the question arises of the structure 
of the superconducting state which not only will have to 
involve the phase coherence of the charge but also of the 
lattice instabilities driving the charge pairing. This im- 
plies cooperative effects acting between such dynamically 
fluctuating local clusters, which we expect to result in 
fluctuating topological structures of the lattice, where in- 
sulating correlations compete with superconducting ones. 

A field theoretical formulation of the Boson-Fermion 
Model, 36 not taking explicitly into account the lattice 
degrees of freedom, has some similarity to the picture of 
a coarse-grained superconducting state of spatially ran- 
domly distributed bosonic modest However, the main 
difference between this approach and the one based on 
the Boson-Fermion Model is that in the latter the local 
inhomogeneities can fluctuate like in Josephson junction 
arrays with inter-grain tunneling. 

Treating the cuprates within this scenario will as a next 
step require a microscopic treatment of the ligand envi- 
ronment. This means: (i) explicitly taking into account 
the different formal valence states of the Cu-ligand com- 
plex, (ii) incorporating the strong Hubbard U correla- 
tions in our present scheme, resulting in local Zhang-Rice 
singlets, and (iii) study the dynamical charge exchange 
between such cation-ligand complexes and the surround- 
ing metallic matrix beyond the single cluster model (dis- 
cussed here) and beyond any adiabatic approximation. 71 
These and other related questions will be dealt with in 
some future work. 
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